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Abstract 

We consider Laguerre polynomials L^ n \nz) with varying negative parameters 
a n , such that the limit A = — \\m n a n /n exists and belongs to (0,1). For A > 1, 
it is known that the zeros accumulate along an open contour in the complex plane. 
For every A € (0, 1), we describe a one-parameter family of possible limit sets of the 
zeros. Under the condition that the limit r = — lim n — log [dist (a n , Z)] exists, we 
show that the zeros accumulate on T r U [/3i, fo] with (3\ and /?2 only depending on A. 
For r £ [0, oo), T r is a closed loop encircling the origin, which for r = +oo, reduces 
to the origin. This shows a great sensitivity of the zeros to a n 's proximity to the 
integers. 

We use a Riemann-Hilbert formulation for the Laguerre polynomials, together 
with the steepest descent method of Deift and Zhou to obtain asymptotics for the 
polynomials, from which the zero behavior follows. 
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1 Introduction 



The zeros of Laguerre polynomials with negative parameters a have a remarkable 
behavior. For a > — 1, these polynomials are orthogonal with respect to the weight x a e~ x 
on [0, oo), and the zeros are all real and positive. This is no longer true if a < — 1. The 



explicit representation, see |31 



k=0 



n + a\ (—z) 



n 



shows that has a zero of order \a\ at if a is an integer in {— n, —n + 1, 
these values of a we have 



(n + a)\ 



(1.1) 

-1}. For 
(1.2) 



This means that in addition to the |a| zeros at 0, there are n — \a\ zeros in (0, oo). For a 
non-integer a G (— n, —1) it is known that there are exactly n — [—a] positive zeros and 
for a < —n there are no positive zeros at all, see |JT], §6.73]. In all cases there is at most 
one negative real zero. 

Non-trivial limiting behavior takes place for the zeros of scaled Laguerre polynomials 



nz) with varying negative parameters a n such that the limit 

a, 



lim — - 

n— >oo fl 



A > 



;i.3) 



exists. The asymptotic location of the zeros depends on A as shown in Figure [I] which is 
taken from see also [25|. See [26] for analogous plots of zeros of Jacobi polynomials 
with negative parameters. 



O 



( 



( 



Figure 1: Zeros of Laguerre polynomials L n n (nz) for n = 40 and A = 0.81 (left), 
A = 1.01 (middle), and A = 2 (right). 

From the plots in Figure [I] it seems clear that the zeros accumulate along certain 
contours in the complex plane. For A > 1, the contour is an open arc, while for < A < 1, 
the contour consists of a closed loop together with an interval on the positive real axis. In 
the intermediate case A = 1, the contour is a simple closed contour. 
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The case A > 1 is well-understood. Indeed, Saff and Varga |29| showed that all of the 
zeros accumulate on a well-defined open arc. Saff and Varga stated their results in terms 
of the zeros of the numerator and denominator polynomials of Pade approximants to the 
exponential functions, but these polynomials can be identified as Laguerre polynomials 
with integer parameters less than —n. Uniform asymptotics for the Laguerre polynomials 



Ln (nz) in case A > 1 were obtained more recently in fli~4" |, ||21|| , and |32 



The case < A < 1 is less well-understood. The asymptotic expansion of Dunster 
TM is also valid for this case, but his results on zeros are restricted to the case A > 1. 



Martinez et al. []25| conjectured that if (a n ) is a sequence of non-integers such that the 



limit ( |1.3|) exists with < A < 1, then the zeros of the polynomials L^ n \nz) accumulate 
on a well-defined contour which they identify as a trajectory of a quadratic differential. As 
stated this conjecture is not true. Indeed, if every a n is an integer, then there are \a n \ zeros 
at the origin. Since the zeros depend continuously on the parameter, it will follow that 
also for a n very close to the integers, some of the zeros of the polynomials L^ n \nz) are 
still close to the origin. However, as we will show, the conjecture in |25[ is true provided 



the distance of the parameters a n to the integers is not exponentially small as n — > oo. 

It is the purpose of this paper to study the possible limiting behaviors of the zeros of 
L [ n n \nz) as n -»• oo such that (|L3D holds with A G (0, 1). We define 



ft = 2 - A - 2y/l-A, (3 2 = 2-A + 2y/l-A. (1.4) 

and 

R{z) ={z- Pi) 1/2 (z - (3 2 ) l '\ zeC\ [/3i,/3 2 ], (1-5) 

where the branches of the square roots are chosen so that R{z) ~ z as z — > oo. We also 
put 

<K*) = l['—d*> * e C \ ((-oo, 0] U [ft, oo)) , (1.6) 

where the path of integration from f3\ to z lies entirely in C \ ((— oo, 0] U [/3i, oo)) except 
for the initial point (3\. 

Due to the pole of R(s)/s at s = with residue R(0) = —A, the function has a jump 
on the negative real axis, 

<f>+(x) — 4>-{x) = — Am, x G (— oo,0), (1.7) 

where <p + and 0_ denote the limits from above and below, respectively. In view of ( |1.7| ), 
Re is a well-defined harmonic function on C \ ({0} U [j3i, oo)). 

Theorem 1.1 For every r > 0, the following hold. 

(a) There exists a simple closed curve T r C C \ ({0} U (fli, oo)) encircling once in the 
clockwise direction, so that 

Re0(z) = r/2 for z G T r . 
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(b) The (a priori complex) measure 



duJs) = as on L r 

v ' 2m s 

where T r is oriented clockwise, is a positive measure with total mass A. 
(c) The measure 



d\x r = dv r + 



2ttx 



-X[Pl,P2]( X ) dx 



is a probability measure on T r U [/?i,/3 2 ]. 
Some of the contours T r are shown in Figure |2| together with the interval fy] 




Figure 2: The interval P2] and some of the contours T r for the value A = 0.99. 

Now we consider a sequence L^ n \nz) of rescaled Laguerre polynomials with varying 
parameters a n such that ( |1.3| ) holds with A 6 (0, 1). We say that the zeros accumulate on a 
compact set K if for every open set Q D K, there is no such that for every n > uq, all zeros 
of l}n n \nz) are contained in Q. We say that a probability measure u is the asymptotic 
zero distribution if 



1 n 

lim - V o 2j = u 



in the sense of weak* convergence of measures, see e.g. |28[| . Here zi 5 „, . . . ,z n>n are the 
zeros of Ln (nz), counted according to multiplicity, and 5 Z is the unit Dirac measure at 
z. 

Theorem 1.2 Let (a n ) be a sequence such that the limit ( |1.3| ) holds with A £ (0,1). 
Suppose that the limit 



L = lim [dist (a n , Z)] 

n— *oo 



1/n 



;i.9) 



exists. Then the following hold. 
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(a) If L = 0, then the zeros of L^ n \nz) accumulate on {0} U [fa, fa] and 

A5 + ^ X[piMK x ) dx I 1 - 10 ) 

is the asymptotic zero distribution. 

(b) If L = e~ r with < r < oo, then the zeros of Ln (nz) accumulate on T r U \fli, fa] 
and the measure \x r given in (|1.8|) is the asymptotic zero distribution. 



Remark 1.3 The case L = 1 in Theorem 1.2 represents the typical case, in the sense that 
if a sequence (a n ) would be chosen randomly according to some reasonable distribution, 
then then 



lim [dist (a n , Z)] 



l/n 



1 



with probability one. So in the typical case, the zeros cluster on [f3i, fa] and along the 
outer curve IV The case of limit L < 1 in ( p_.9[) is more special, since for that to happen 
the parameters a n should be very close to integers. 

To illustrate this, we have plotted in Figure |3] the zeros of Ln _nA ' ) with n = 40, for the 
values A = 0.81 and A = 0.799999975. Then the distance of a n = —nA to the integers 
is 0.4 and 10 -6 , respectively. There is a clear separation in the plot on the right between 
the zeros on the interval \j3i, fa] and the cluster of zeros around the origin. For the plot on 
the right, fa = 0.31 • • • and fa = 2.09 • • • . The rightmost zero in the cluster has real part 
0.14- ••. 



Figure 3: Zeros of Laguerre polynomials Ln '(nz) for n = 40 and A = 0.81 (left), 
A = 0.799999975 (right). The two cases correspond to dist (a n , Z) = 0.4 and dist (a n , Z) = 
10~ 6 , respectively. 



Remark 1.4 Theorem 1.2 also applies to subsequences. So, if (a n ) is a sequence such 
that the limit (|1.9| ) does not exist, then we can pass to a subsequence such that J1.9Q 
exists along that subsequence. Then we find an asymptotic distribution of zeros for any 
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subsequence such that the limit ( p..9|) exists. By constructing special sequences we can 
arrange it that different subsequences have different asymptotic zero distributions, and it 
is even possible to construct examples in which every measure \x r appears as the asymptotic 
zero distribution along some subsequence. 



Remark 1.5 For fixed n and decreasing a < 0, the dynamics of the zeros of L„ is similar 
to the dynamics of the zeros of certain hypergeometric polynomials studied recently by 
Driver and Duren [12 , |i~3f . For the Laguerre polynomials, the zeros are all real and positive 



if a > — 1. For a = — 1, there is one zero at the origin. This zero becomes negative if a 
becomes smaller than —1. As a decreases, the smallest zero reaches a minimum and then 
starts to increase again. For a = —2, it is at the origin again, and simultaneously the 
second smallest zero has reached the origin. A collision takes place, and for a somewhat 
smaller than —2, the two zeros separate and leave the real axis, moving into the complex 
plane. When a has reached the next integer value —3 they both return to the origin, and 
a third zero coming from the positive real axis joins them in a three-fold collision. As a 
decreases from —3 to —4, the three zeros move into the complex plane, one on the negative 
real axis, and two away from the axis, coming again to the origin for a new collision at 
a = —4, where they are joined by another zero coming from the positive real axis. And 
so on, until all n zeros have left the positive real axis. At a = —n, all n zeros meet at the 
origin. For a < —n, they separate again, and they all move into the complex plane. If n 
is odd, then one zero remains on the negative real axis. 



Remark 1.6 In statistical literature, see IU] and the references cited therein, the mea- 
sure (|1.10| ) is known as the Marchenko-Pastur distribution [24|. It appears as the limiting 



spectral distribution of Wishart matrices, which are random covariance matrices ^XX T 
where X is an n x p matrix with independent, identically distributed entries. The measure 
( P-.10I ) corresponds to the limit n,p — > oo, n/p — > 1 — A e (0, 1). 



In the theory of free probability, the measure ( |1.10 ) is also known as the free Poisson 
PI 



distribution, see 



and the references cited therein. 



We prove Theorem 1.1 in section 2. The rest of the paper is devoted to the proof of 
Theorem 1.2. It is based on a characterization of Laguerre polynomials by a Riemann- 
Hilbert problem [IS], |2"I|, and an asymptotic analysis of the Riemann-Hilbert problem with 
the steepest descent / stationary phase method introduced by Deift and Zhou 0, and 
developed further in |4|, || [7|, || |20|, p2]| . See |5j for an introduction to this method, with 
a special emphasis on the applications to orthogonal polynomials and random matrices. 
In [^TJ we used the Riemann-Hilbert method to study Laguerre polynomials with varying 
negative parameters a n such that lim n a n /n < —1, which is the case A > 1. In this case, 
the zeros accumulate along an open contour as shown in Figure 1. The determination of 
this contour is a first step in the analysis in ||21|| . Riemann-Hilbert analysis in which a 
contour selection method was required appeared also in the recent papers |l|, 



19 . 
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In the Riemann-Hilbert method, the endpoints of the contour need special attention. 
Following H |7j , the local analysis around the endpoints is done in |21] with the help of Airy 



functions by constructing a so-called Airy parametrix. The analysis for the case < A < 1 
proceeds along similar lines. However, the fact that the zeros accumulate along a more 
complicated contour causes some new features. In particular we will be dealing with a 
closed loop r together with an interval [fa, fa] on the real line. The point fa where the 
closed loop and the interval intersect is a new special point. The surprising fact, however, 
is that after suitable transformations, the transformed Riemann-Hilbert problem is very 
similar to the transformed Riemann-Hilbert considered in [ff| . It then follows that the local 
analysis around fa can be done with the Airy parametrix as well. 

The Riemann-Hilbert method gives strong and uniform asymptotics of the Laguerre 
polynomials Ln' n \nz) in the whole complex plane. This is really much more than what 
is needed for the proof of Theorem 1.2. Since we do not have an application for strong 
asymptotics (although we could envision there is one), we have not pursued this matter 
any further. As a sample of the kind of results that could be obtained, we state here the 
asymptotic result in the oscillatory region /? 2 ) on the real line. Suppose that the limit 
( |OD holds with A 6 (0, 1). For every n G N, we put A n = so that A n — > A as 

n — > oo. Let fa >n and fa ;U be the values (|1.4j) corresponding to A n instead of A. Then we 
have, uniformly for x in compact subsets of (Pi, fa), 

L^\nx)- 2 ' 



77! 



COS 777T 



fas 



\J fan 


~ fan (fan ~ X 




fa,n)(fa,n ~ S) 


2tts 




~ fa,n)(fa,n - s) 


2TXS 



1 . 2z-(fa >n + fa >n ) 
- arcsm 

2 fa,n — fa.n 




: sh, ( nn / ^ ' ds - \ arcsin 2 " ± ^ I O I - ] I I II ! 



2 Proof of Theorem 1.1 

Proof. The curves where Re <p(z) is constant are known as trajectories of the quadratic 
differential 

_m#^ = _(.- ft )(,-A )<faa 

which has simple zeros at fa and fa and a double pole at 0, see |30|]. From the local 



structure at a simple zero of the quadratic differential, see |3(], §7], we know that three 
trajectories emanate from fa at equal angles 2n/3. One of these is the interval [fa, fa]. Let 
T be the trajectory emanating from fa at the angle — 2tt/3, and let z(t) be the arclength 
parametrization of r with z(0) = fa. So Re [0(z(i))] = for every t, which implies in 
particular that T does not approach the origin. Since r leaves fa at an angle — 27r/3, 
we have for small t > (and thus for z(t) close to fa), that \z(t)\ decreases as t increases. 
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Then \z(t)\ continues to decrease until for some to, we have 



Re 



z'jto) 



(2.2) 



Since Re [(j)(z(t))] = for every t, we find by differentiating ( |f .6| 



Re 



z'(t) 



Rim 

z(t) J 







(2.3) 



for every t. Thus for t = to we have from ( |2.2j ) and 
is only real for z real, it follows that z(to) is real, So 
axis. Since z(0) = Pi, we have |^(to)| < Pi- From ( 
strictly positive for x G (0, Pi), since R(s) < for s 
have a curve in the lower half-plane from Pi to z(to 
also have Re = on the mirror image of this curve, 
form the full trajectory Tq which encircles the origin, 
for r = 0. 

See Figure f| for a plot of T for two values of A. 
where Re0(^) = (solid lines). The dotted lines are 



(gj) that R(z(t )) e E. Since R{z) 
\z(t)\ decreases until To hits the real 
|1.6| ) it is clear that <ft(x) is real and 
G (— oa,Pi). Thus z(to) < 0. So we 
) where Re0 = 0. By symmetry, we 
Then the curve and its mirror image 
This proves part (a) of Theorem 1.1 

Also plotted are the other contours 
the curves where either d>(z) or 



are real. 



R(s) 



ds 



(2.4) 






Figure 4: Curves where <f) is purely imaginary (solid lines), or where either or are real 
(dotted lines). The closed loop is T . The values for A are 0.99 (right) and 0.81 (left). 



Let flo be the region that is encircled by r . Then Re0 is harmonic on Q \ {0} with 
zero boundary values on r . From (|1.6|) it follows that Re <p behaves like ^ log \ z\ as z — > 0. 
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Thus (2/A)Re<j) is the Green function of Qq with pole at 0. Since To is a simple closed 
curve, it then follows that for every r > the level curve 



r r = {z G Q | Re0(Y) = r/2} 

is a simple closed curve encircling the origin. This proves part (a) for every r > 0. 
For part (b) we first note that by Cauchy's formula 

duJs) = — ! ^ds = -MO). 
2ni J Vr s 

We have a —sign since r r is oriented in clockwise direction. Since R(0) = —A, we get 
J r dv r = A > 0. If z(t) is arclength parametrization of T r with z(0) = fii, then 

Mt) 1 ,*(t) R , s j 

1 1 — ^ / az/ r = / as 

7/3i 27rz J Pi s 

is differentiable with derivative ^- ^^ ^(t), which is real by the definition of T r , and 

non-zero since R is non-zero on T r . So dv r can only increase along T r and it follows 
that v r is a positive measure with total mass A. This proves part (b). 

Then for part (c) it is enough to observe that 



which can be shown by direct calculation, or by a contour deformation argument. □ 

3 Proof of Theorem 1.2 
3.1 Preliminaries 

We assume that the sequence (a n ) is as in Theorem 1.2. We put 

4. = -—- (3-1) 

n 



so that lim A n = A. Since A e (0,1), we may assume without loss of generality that 

n— *oo 

A n G (0, 1) for all n, so that a n G (— n, 0). We may also assume that a n ^ Z. Indeed, if 
a n G Z, then by (|1.2j) the zeros of Ln (nz) are the zeros of L^~+"^(nz), together with a zero 
at of multiplicity \a n \. From the asymptotic behavior of zeros of Laguerre polynomials 
with varying positive parameters, see [11], |1^, [17] ^3|, part (a) of Theorem 1.2 follows, in 



case the a n are integers. So we may and do assume that a n ^ Z. 
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For every n, we use a subscript n to denote the notions ( |1.4j) , (|1.5| ), ( |1.6| ), but with A 
replaced by A n . Thus 



(3 1)n = 2-A n -2^Jl-A n) fa >n = 2 - A n + 2^1 - A„. 
/2 n (z) = (s - /3 lin ) 1/2 (z - &, n ) 1/2 , 2 e C \ [A, n , /3 2 , n ], 



and 



We also put 



<Pn{z) 



/3i, 



0n(^) 



1 



2 G C \ ((-oo, 0] U [p 1>n , oo)) . 



ds, z G C \ (-oo, f3 2 , n ] 



(3.2) 
(3.3) 

(3.4) 
(3.5) 



where the path of integration from p 2 ,n to z lies entirely in C \ (— oo,/3 2 , n ], except for 
the initial point ^ 2 ,n- We denote by Iq the simple closed contour around given by 
Re<j) n (z) = 0, whose existence is guaranteed by Theorem 1.1(a). 



3.2 The Riemann-Hilbert problem 

For every unbounded contour £ in C \ [0, oo) from a point +oo — iy to +oo + iy for some 
y > 0, we have the orthogonality relation 

J L { - nA "\z)z k z- nA "e- z dz = 0, for k = 0, 1, . . . , n - 1, (3.6) 

where the branch of z _nA ™ is taken with a cut along the positive real axis. In addition, 
since nA n ^ Z, 

L^ nAn \z)z n z- nAn e- z dz ^ 0. (3.7) 

The relations (|3.6|) and ( |3.7| ) were proved in Lemma 2.1]. They follow from the 
Rodrigues formula for Laguerre polynomials. 
The rescaled monic polynomial 

Pn{z) = j^I^H (3.8) 

is then characterized by a Riemann-Hilbert problem on a contour £ as shown in Figure [5], 
oriented from +oo — iy to +oo + iy. For any oriented contour £, we follow the usual 
convention that the +- side is on the left and the — side is on the right while traversing the 
contour. We say that a function Y defined on C \ £ has boundary values on S if for every 
zGS, the two limits 



Y+(z) = lim Y(z'), Yl(z) = lim Y(z') 

+side9z'— *z — sidc3z'^z 



exist. 
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Figure 5: Example of contour E 



Riemann-Hilbert problem for Y: 

Let E be as above. The problem is to determine a 2 x 2 matrix valued function Y : C \ E — 
C 2x2 such that the following hold. 

(a) Y(z) is analytic for z G C \ E, 

(b) Y(z) possesses continuous boundary values for z G E, denoted by Y + (z) and Y-(z) 
and 

(i 7 — nAn p —nz 
Q x )■ for.:GE. 

(c) y(z) has the following behavior as z — > oo: 
y(*) = (j + O (T^ ( ^ ) a .s : oc. : G C \ E 



The Riemann-Hilbert problem has a unique solution, see [f21~1 , Proposition 2.2], given by 

(3.9) 



Y z) 



C-* 

\ 2J / 

where (J n _i(z) = d n L^l"' A ™'' (nz) is a polynomial of degree n — 1, and the constant d n is 
chosen such that 122(2) behaves exactly like z~ n as z — > 00. 

Anticipating the distribution of zeros as shown in Figure [I] for the case < A < 1, we 
modify the Riemann-Hilbert problem for Y as follows. We let part of the contour E come 
to the positive real axis and put 

S = E^ = r{ ) n) U[/3 1 , n ,oo). 
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The curve is oriented in the clockwise direction, and [fli >n , oo) is oriented from left to 
right. The unbounded component of C\S^ n ) is denoted by and the bounded component 
by Q as indicated in Figure ||. 




Figure 6: Contour E^ = U \p 1>n , oo) 

Then it immediately follows that P n is characterized by the following Riemann-Hilbert 
problem. 

Modified Riemann-Hilbert problem for Y: 

Let E 1 -™) be as described above. The problem is to determine a 2 x 2 matrix valued function 
Y : C \ E -> C 2x2 such that the following hold. 

(a) Y(z) is analytic for z E C \ 

(b) Y(z) possesses continuous boundary values for z E E™, denoted by Y + (z) and Y-(z), 



and = Y^(z)Vy(z) where 

1 



for zerf 1 , 



(3.10) 



Vy(z) 



1 c n z 




and c n = 2i sm(nA n 7r) . 



(c) Y(z) =1 + 



z n 
z~ r 



for z E (/3i n ,oo), 



asz-> oo. 



(3.11) 



Then as before we have that the modified Riemann-Hilbert problem for Y has a unique 
solution and that 



Y 11 (z)=P n {z). 



(3.12) 
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The constant 

c n = 2i sm(nA n n) (3.13) 

appearing in part (b) will continue to play a role later on. Since a n is not an integer, c n is 
different from zero, and it can be used as a measure for the distance of a n to the integers. 
In fact, we have that 

L = lim |c n | 1/n , (3.14) 

n^oo 

where L is given by ( |1.9|) . 



3.3 First transformation Y \—> U 

The first transformation of the Riemann-Hilbert problem is based on the construction of 
a so-called g-function. 

Definition 3.1 We define the function g n by 

9n{z) = ^-l \og(z-s)^^-ds, zeC\SW. (3.15) 

For s G [/3i, n , /52,n], the logarithm log(z — s) is taken with a cut along [s, oo), for s G 
the branch of the logarithm log(z — s) is chosen with a cut along the part of Iq that starts 
at s continues along Iq to the right until /3 1)n , and then consists of [/3i,n, oo). 

Thus g„ is the log-transform of the probability measure [j,q, see (P-.8|), but corresponding 
to the value A n instead of A. Note that g n is defined and analytic in C \ S^ n \ 

The following proposition describes the jump properties of g n . They follow by taking 
into account the various branches of the logarithm, see also [^TJ. We use C + and C~ to 
denote the upper and lower half-plane, respectively. 

Proposition 3.2 (a) There is a constant i n such that 

1 ( ^\A n iri-(j) n {z), for z G ^nC* 

g n (z) = -{A n logz + z+£ n ) + < (3.16) 

2 { ±±A n m + <j) n (z), for z G ^ H C±. 

Here log z = log \z\ + i arg z with arg z G [0, 2ir) . The number £ n is explicitly given by 
£ n = 2g n (x n ) — (A n \ogx n + x n ), where x n is the intersection o/Tq with the negative 
real axis. 

(b) We have the following jump for g n on , 

-2cj) n (z) =F A n m for z G T^ l) n C ± , 

-2<p nt+ (z) - 2A n 7ii = 20 n _(z) - 2A n vri for z G (/?i, n , /3 2>n ], 

-2ni for z G [/3 2 ,n, oo). 

(3.17) 



9n,+( z ) - 9nA z ) 
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(c) We have, with the same constant £ n as in part (a) 

A n logz + z + £ n , 

A n log z + z + A n iii + £ n , 

k A n log z + z + A n iri + £ n - 



9n,+(z) +g n -{z) 



Proof. Taking the derivative of ( 3.15|) we get that 

; iv 



9n( Z ) 



Itxi /pC») z — s s 2m 





1 R r , 



£ ^ -n( n ) 

for z £ Tq , 

/or Z G [Pl,n,p2,n\, 

for z G (/3 2 ,„, oo). 

(3.18) 



-ds. 



\PiM 



The two integrals can be calculated using contour deformation and Cauchy's formula. The 
result is that 



9n( Z ) 



Rn(z) |_ An _[_ I 



1 f Rn(z) | A n | ^ ^ 



for Z G f^oo, 

for z G ^o- 



(3.19) 



We integrate g' n from x ra to z along a path in C \ (£( n ) U R), except for the initial point x n . 
We obtain from ( |3.19| ) and the definition of £ n , 



9n(z) = g n (x n ) + / g' n (s)ds 



-(A n \ogz + z + 




1 r* Rn(s) 



ds for z G Vt c 



_1 T2 

2 Jx 



+ 1 f* f or 2 e ft 

2 J x n s u ' 



where we also used the fact that g„ >+ (x„) = g n -{x n ) because of the choice of the branches 
of the logarithm in ( |3.15| ). [Here g n>+ (x n ) and g n -(x n ) denote the limiting values of g n (z') 
as z' — > x n from within and Q , respectively.] By the definition (|3.4| ) of 4> n (z), we have 

— — -(is = 0(2) — <fi±(x n ) for z G C 1 * 1 , 



2 s 



and ( p.lfip follows, since <f>±(x n ) = ^^A n 7ci. This proves part (a) of the proposition. Parts 
(b) and (c) follow immediately from part (a). □ 

With the o-function, the constant £ n , and the number c n = 2ism(nA n ir), we perform 
the first transformation of our Riemann-Hilbert problem. We define for z G C \ T,( n \ 



U(z) = (c n )"5 CT 3 e -|™^a3y^) e -n« ? „Wa3 e |n^ ff 3( Cn )|a3 > 



(3.20) 



Here 03 is the Pauli matrix 03 



'1 

From the Riemann-Hilbert problem for Y and the jump relations satisfied by g n , it 
follows by a straightforward calculation that U is the unique solution of the following 
Riemann-Hilbert problem. 
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Riemann-Hilbert problem for U : 

The problem is to determine a 2 x 2 matrix valued function U : C \ E^ — > C 2x2 such that 

(a) U(z) is analytic for z E C \ E^, 

(b) U (z) possesses continuous boundary values for z E S^ n \ denoted by U+(z) and U-(z), 
and U+(z) = U-(z)Vu(z) where 

( 2ncf>„ t+ (z)+2nA 
e 2n^.(z)-2nA n ni ) for Z E (/?l,n, 02,n]i (3-21) 

/ 1 „-2n0 n (z) \ 

Vu(z)=l L Q e 1 J for z E (^2,n> oo); (3.22) 

/ „2n0„(2)±n J 4 n 7ri _— 1 \ 

= n p -2n M W- for z G fl C ± . (3.23) 



(c) U(z) =I + oO^ as z -> oo. 

3.4 Second transformation U *—>T 
We choose a small e > 0, so that the two disks 

A £ 0%) = {z G C | \z-/3j\<e}, j = 1,2, (3.24) 

are disjoint. Then there exist no such that for every n > uq and j = 1,2, we have that 
Pj jH E A £ ([3j) and that the contour E^ is within an e-neighborhood of E = r U [/3i , /?2] ■ 

Now for n > no, we open up regions around r o n) and [f3[ n) , (3^}. The regions are called 
Jl^, Qf, ^2, and as shown in Figure 0. The exterior region is called f2^. The boundaries 
of the regions form the contour [xq, oo) U U^=o ^j-> wnere is the intersection of Si with 
the negative real axis. The contour is oriented as in Figure [?]. We can (and do) take the 
contours Ej to be independent of n outside of the disks A e (/3j-) and so that Re (ft < on 
the contours E^ outside the disks. We can clearly arrange for that. Indeed the solid lines 
in Figure |] are the curves where Re0 = 0, and this includes the closed contour r . The 
region where Re < is the unbounded white region shown in Figure []. The dark region 
is the region where Re0 > 0. So away from the points (3\ and fa, we choose the contours 
Ej in the white region. 

We also take E 3 - for j = 1, . . . , 5 to be within distance < 2e from E , and such that the 
point where Ei, E 2 , and E 4 meet, and the point where Ei, E 3 , and E 5 meet are outside 
the disks A e (/?_,■)• 
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-0.5 0.5 1 1.5 

Figure 7: Contour for the Riemann-Hilbert problem for T, for the value A n = 0.99. The 
contour consists of a semi-infinite real interval [xq,oo), the closed curve Tq, and the five 
arcs Ej, j — 1, . . . , 5. The contour bounds the regions VLq , fif , f2 2 , f2 3 , and 



0.8 - 



0.6 - \ \ 

\ \ 




-0.6 - / 

l 

' / 

-0.8- / / / 

/ / 



-0.5 0.5 1 1.5 

Figure 8: The dark region is the region where Re0 > 0, the white region is the region 
where Re0 < 0. The parts of the contours for j = 1, . . . , 5 that are outside the disks 
are chosen in the white region and within distance 2e from Eo- The particular plot is for 
the value A = 0.99. 
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With this choice of contours we define T for n > no, by 

T(z) = U(z) ( ^ e . 2n ^lnA nm ) for z G fi± (3.25) 

T(z) = U(z)( _ Cne 2nMz)±nA n ni J ) for z G fif , (3.26) 

T(z) =E/"(*) ^ _ e 2n0 n J)+2nA n ^ J J for z G ft 2 , (3.27) 

= U(z) ( e2n4>n {z l 2nAnni J) for^Gfis, (3.28) 

T(z) = for 2 G (3.29) 

In ( |3.26|) we have defined T with a different formula in Qf and fijf. So one might expect 
a jump across fl± D R. However, we have as in (|1.7| ) 

n ,+ (x) = <j) n -(x) - A n ni, for x < 0, (3.30) 



so that ( ggp gives that T+(a;) = T_(x) for a; G fii DR. Similarly, by O) and (|3T30| ), we 



have T" + (x) = Tl(x) for a; G f2o H R . It follows that T has an analytic continuation across 
(x ,0). 

Next, we see by direct calculation from (|3.23|) and (|3.25|) that T + = T_ on and 



therefore T has an analytic continuation across as well. So T is analytic on the 
complement of the reduced contour S T shown in Figure |9|. The regions bounded by the 
contour E T are denoted by Qq, f2 2 , ^3, and fl^ as indicated in Figure || 

Now from the Riemann-Hilbert problem satisfied by U, we obtain by straightforward 
calculations the jumps for T on S T . The result is that T is the unique solution of the 
following Riemann-Hilbert problem. 

Riemann-Hilbert problem for T: 

The problem is to determine a 2 x 2 matrix valued function T : C \ S T — > C 2x2 such that 
the following hold: 

(a) T(z) is analytic for z G C \ £ T , 

(b) T(z) possesses continuous boundary values for z G S T , denoted by T + (z) and T_(z), 
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0.4 - 




-0.8 - 



-1 ' — 
-0.5 



Figure 9: Reduced contour S T for the Riemann-Hilbert problem for T, for the value 
A n = 0.99. The contour consists of the semi axis [0, oo) and the five arcs Sj, j — 1, . . . , 5. 
The contour bounds the regions Qq, f^, ^3, and fl^. 



such that T + (z) = T-(z)Vr(z) where 
V T (z) = 



V T (z) = 
V T {z) = 
V T (z) = 
V T (z) = 
V T (z) = 
V T (z) = 



1 

c e 2n<f> n (z)±nA n m j 



1 

p 2ntj> n {z)+2nA 



1 







2n<t> n (z)— 2nA n ni -y 



1 

e 2n<t> n (z) j 
^ e -2ra</> n (z) 

1 



1 
-1 



^ e -2n4> n (z) 

1 



for zeSifl C ± , 

for z G S 2 , 

for z G S3, 

for z G £ 4 U £5, 

for z G [0,/3i,„), 

for z G (P l7 n,P2, n ), 

for 2; G (#2,n, 00). 



(3.31) 
(3.32) 
(3.33) 
(3.34) 
(3.35) 
(3.36) 
(3.37) 
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T(z) = I + O I - ) as 2 oo. 



Remark 3.3 The contour E has been constructed so that Sj \ (A £ (/? 1 ) U A £ ((3 2 )) is n- 
independent, and so that 4>(z) < there. Since <f) n (z) — > (f)(z) as n — > oo, we see that the 
jump matrices in ( |3.31| )-( j3.34| ) converge to the identity matrix as n — > oo. The convergence 



is uniform for z outside the disks A e (/3j). 

Similarly, since <f>(z) > for z G [0, > for 2: G (/3 2 ,oo), and 0„ — > (f>, and 

</>„ — > (f) as n — > 00, the jump matrices in (|3.35|) and ( p . 3 7|) also converge to the identity 



matrix as n — > 00, and the convergence is uniform for z outside the disks A £ ((3j). 
3.5 Construction of the parametrix for T 



Observe that we may write the jump matrices in ( [3.321 ) and ( 3.33|) also as 



1 

e 2n4> n (z) y 



Then it follows that we have a situation which is basically the same as the one considered 
by Deift in 0, p. 199]. There the part Si is absent, and the contour extends to —00 along 
the negative real axis. However, these differences are not important since on these parts 
the jump matrices are uniformly / + 0(e~ cn ' 2 ') asn-> 00, for some c > 0. 

It follows that we can do the same analysis as in || (see also @, 0) to construct the 
parametrix for T . The behavior away from (3\ and (3 2 is governed by the solution of the 
following Riemann-Hilbert problem. 

Riemann-Hilbert problem for N: 

The problem is to determine N : C \ \Pi >n , /?2,n] ~~ > C 2x2 such that the following hold. 

(a) N(z) is analytic for z G C \ [/3i )Tl , 02,n\, 

(b) N(z) possesses continuous boundary values for z G (Pi, n , P2,n) , denoted by N + (z) 
and N_(z), such that 

N+(z) = N-{z)( ° M forzG(/3 1)n ,/3 2i „), (3.38) 



•1 t 

and N has at most 4-root singularities at /3i jn and p 2 ,n, 

(c) N(z) =1 + 0(1) as z °°- 

Writing a„(z) = [ z ~^ 2 ' n l/, d , we then have that the solution of the Riemann-Hilbert 
problem for N is 

q n (2)+q n (2:)~ 1 a, 1 (z)-q n (2)~ 1 

JV(«)=| I for 2 GC\D9 1>n , /%,„]. (3.39) 

-2i 2 
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Near 0% and /3 2 we construct local parametrices, denoted by P\ and P 2 , respectively. 
Recall that we have chosen e so that the disks A e (/3i) and A £ (/3 2 ) are disjoint. 



Riemann-Hilbert problem for Pj, j = 1,2: 

The problem is to determine matrix valued functions Pj : A £ ((3j) \ T, T — > C 2x2 such that 
forj = 1,2, 



(a) P/(z) is analytic for z G A £ (/3j) \ X 3 

(b) Pj satisfies the same jump conditions on S T fl A £ (/3j) as T does. 

(c) We have as n — > 00, 

= LT + O (^j J N(z) uniformly for \z - fa\ = e. (3.40) 

The Riemann-Hilbert problems for Pi and P 2 are solved as in [||, |7| with the use of Airy 
functions, see also || [21]. Note that Pi and P 2 , as well as N, depend on n. 



The formula for P 2 is as in [0, Eq. (7.24)]. It uses a conformal mapping 

C = /n(z) = ( -Z n <t>n{z) 



from A £ (/3 2 ) onto a neighborhood of in the £-plane, and an explicit 2x2 matrix ^(C) 
built out of Airy functions which is given in p. 1522]. The superscript a denotes an 
angle that may be any number in (tt/3, it), possibly be n-dependent. Then P 2 has the form 



P 2 (z) = E n {z)W{f n {z)) 



e n<t> n (z) Q 

Q e -n-4>n(z) 



where E n is an analytic pre-factor that takes care of the matching condition Q3.40| ), see |7|, 



Eq. (7.23)] for the precise form of E n . The construction of the parametrix Pi near the left 
end-point fix is similar, see [[7|, pp. 1526-1527]. 

3.6 Final transformation T 1— > S 

Having N, Pi, and P 2 , we define 



S(z) = T(z)N(z)- 1 for^eC\(E T UA7p 1 )UA £ C^)), (3.41) 
S(z) = T(z)P 1 (z)- 1 hrze A £ (A)\S T , (3.42) 
S(z) = T(z)P 2 (z)- 1 for z G A £ (/3 2 ) \ S T . (3.43) 
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0.8 
0.6 



-0.2 
-0.4 
-0.6 
-0.8 




-0.5 



0.5 



1.5 



Figure 10: Contour S 5 for the Riemann-Hilbert problem for S, for the value A = 0.99. 
The two circles are the boundaries of two ^-neighborhoods of (3\ and j3 2 , with e = 0.1. 

Then S is defined and analytic on C \ (S T U dA £ (/3i) U dA £ ([3 2 )). However it follows from 
the construction that S has the identity jump on the interval (f3\ + e, fa — e) and on the 
parts of S T that are within the disks A e (/3j). Therefore S has an analytic continuation to 



C \ X , where S is the system of contours shown in Figure 0. The contour £ does not 
depend on n. 

Then S satisfies the following Riemann-Hilbert problem. 
Riemann-Hilbert problem for S: 

The problem is to determine S : C \ S s — > C 2x2 such that the following hold. 

(a) S(z) is analytic for z G C \ S 5 , 

(b) S(z) possesses continuous boundary values for z G S 5 , denoted by S+(z) and ^-(z), 
such that S+(z) = S-(z)Vs(z) where 



V s (z) = P 1 {z)N{z)- 1 
V s {z) = P 2 (z)N(z)- 1 

V s (z) = N^Vr^Niz)- 1 



for dUe(pi), 
for dU £ (/3 2 ), 

foTzeE s \(dU E (fa)UdU £ (p 2 



(3.44) 
(3.45) 

(3.46) 



(c) S{z)=I + 



as z 



oo. 
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All jump matrices for S are close to the identity matrix if n is large. The jump 
matrices in ( p.44 ) and ( |3.45|) are l+0(l/n), and the jump matrices in (|3.46|) are uniformly 
I + 0(e~ cn ) for some c > 0. In addition, the jump matrix ( p. 46 ) on the unbounded interval 
(@2 + e i 00 ) is uniformly / + 0(e _cn ' 2 ') for some c > 0, so that the jump matrix is close to 
the identity matrix both in L2-norm and in Loo-norm on the contour E s . Also note that 
by construction the contour E does not depend on n. Then it follows as in || that 



S(z) 



Oil 

n 



uniformly for z 6 C \ E' 



(3.47) 



as n 



00. 



3.7 Proof of Theorem 1.2 

Tracing back the steps Y \— > U 1— > T 1— > S and using ( |3.47| ) we obtain strong asymptotics 
for Y in every region in the complex plane. In particular since 

Y n (z)=P n (z) = -^L<r^\nz), 

this yields strong asymptotics for the generalized Laguerre polynomials Ln n {nz). 

The formulas we obtain in various regions are similar to the ones obtained in for 
orthogonal polynomials with respect to exponential weights on R. We will not give all these 
results here, but restrict ourselves to the ones that are needed for the proof of Theorem 
1.2. 

Proof of Theorem 1.2. As noted before, we may and do assume that a n (jL Z for every 
n. For every n, we can then set up the Riemann-Hilbert problem for Y, and by ( p.9|) and 
( Pq ) we have 

P n (z) =Yu(z) =Un(z)e n9 " (z) . (3.48) 

Using the formulas ( |3.25| )-( |3.29| ) which express T in terms of U, we find P n in terms of the 
entries of T. The precise formula depends on the region. 

For z in the outer region f^, we find by ( |3.29| ) and ( [3.48] ), 

P n {z) =T ll (z)e n9n{z) . (3.49) 

If, in addition, \z-(3 2 \> e, then T = SW by ( ggTfl and S = I + 0(l/n) by OOTQ , so that 



P n (z)e-^^ = ^(z) ^1 + 



(3.50) 



uniformly for z G \ Ag^). We have written jV( n ) to emphasize the n-dependence of 
iV. As n — >• 00, we have by (|3.39|) 

hm N$\z) = hm a ^ +a ^~ 1 = <M±^T1 (3.51) 
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uniformly for z G \ A e (/32), where a(z) = |^_^]i/4 ■ Since the right-hand side of fl3.51| ) 



does not have zeros in C, it follows from (|3.50 ), ( 3.51 ), and Hurwitz' theorem, that there 
are no zeros of P n in \ A £ (/? 2 ), if w is sufficiently large. From a contour deformation 
argument as in 0, p. 1532], we get that fl3.50| ) holds for z G (/3 2 + £, oo) as well. Thus, by 
the same argument, there are no zeros in this interva l if n is sufficiently large, and so, for 
large n, all zeros of P n are in CIq U O2 U ^3 U A e (/3 2 ). 

Next, we have by fl33Bp and (ft^gQ 



p n (z) = f/ n ( 2 )e n9 " (z) = [ e - 2n *"W :Fni4n7rt rii(z) + c n T 12 (^)] e n9n{z) for 2 G flj, (3.52) 



and by (p|) and (ggp 



P n (z) = f/ n (z) e n9 " (z) = [T n (z) + c n e 2n ^ (z)±nA "^T 12 (2)] e n9 " (z) for z G ftf. (3.53) 



For 1 2; - ^1 > e, we have T = SW by ( ggg ) and S 1 = 7 + 0(l/n) by (|3T47D . Then we get 
from (|3321 ) 



+ 



(«), 



oil 

n 



n 



(3.54) 



uniformly for z G fio \ A e (/?i), and from (|3.53|) 



P n (z)e 



-ng-n(z) 



N$\z) + c n e 2nMz)±nA " m N^\z) 



1+0 
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+ 



7? 



(3.55) 



uniformly for z G fli \ A £ (/5i). 
As n — > 00, we clearly have 

lim iV (n) ( 



a(z)+a(z) 1 a(z)— a(z) 1 
2 2l 

a(z)— a(z) _1 a(,z)+a(z) — 1 



-2i 



from which it follows that the entries of are uniformly bounded and uniformly bounded 
away from zero on any compact set which does not contain f3\ and /3 2 . We also recall that 



L = lim \c n \ 1/n . 

n^oo 

Now we distinguish the cases L = and L = e~ r as in Theorem 1.2. 



(3.56) 
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(a) Case L = 0. 

Let K = Qq \ (A E (fa) U A £ (0)). Then the functions (fi n are uniformly bounded on 
K, so that we get from ( |3.56|) and L = that 

lim | Cne 2n0n(^)TA„«| 1 /" = o un if orm l y f or z e K. 

Since the entries of are uniformly bounded and uniformly bounded away from 
0, it then follows from (|3|) and ( ET55D that 

aJz) + aJzY 1 ( ' . „ /l 



7? 



uniformly for 2 G if fl fi , and 



2 \ \n 

uniformly for z G K r\Q\. Thus for large n, there are no zeros in K, so that all zeros 
are in A e (0) U A £ (/?i) U fi 2 U fi 3 U A £ (fa). Since £ can be chosen arbitrarily close to 
0, and the domains f2 2 and Q3 lie within the 2e neighborhood of [fa, fa], it follows 
that the zeros accumulate only on {0} U [fli, fa]. 

(b) Case L = e~ r with < r < 00. 

To locate zeros of P n we have to balance two contributions in ( |3.54 ) or (|3.55|) , which 



can only be done near the curve T r where Re0 = r/2. Indeed, if K is 

K = {z G T% I dist (z, r r ) > e} 

then there is 5 > such that \<fi n (z)—r/2\ > 5 for z G K and for n large enough. Then 
either c n N^{z) or e -^ n 4>n{z)^nA n m ^ d ornma tes the other term in the asymptotic 
formula ( 3.54Q . The former dominates if Ke<j) n (z) > r/2 + 5 (which is for z in the 



interior region bounded by r r ), and the latter dominates if Re (j) n (z) < r/2 — 5. Then 
it follows that for large n, there are no zeros in K fl Qq. Similarly, it follows from 
( |3.54 ) that no zeros in K fl Q\ for large n. Thus the zeros in Q$ \ A £ (/3i) are in a e 



neighborhood of T r if n is large. In addition, there can be zeros in A £ (fa) U A £ (/? 2 ) U 
f?2 U ^3. Letting e — > 0, we see that zeros can accumulate on T r and on the interval 
[ft, ft] only. 

So now we know that the zeros accumulate on {0} U [fa, fa] in case L = 0, and on 
r r U [fa, fa] in case L = e~ r > 0. It remains to determine the asymptotic zero distribution. 
To that end, we note that we get from ( |3.50| ), 

lim (-\og\P n (z)\-Reg n (z)] = 0, (3.57) 



n— >oo \ fl 
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uniformly for z in compact subsets of f^. From the definition of g n , see ( |3.15| ), it follows 
that 



lim Reg n (z) = / log \z — s\ d[i (s) (3.58) 

n— >oo / 



uniformly for z in compact subsets of C \ (r U [fa, fa]). Here fi is the measure on T U 
[fa, fa] given by Q with r = 0. Thus by (j337|) and flOp , 



lim -log|P n (z)| = / log \z - s\ dfj, (s), (3.59) 

n— >oo j 

for z in f^. 

From the definition of the measures /i r in (|1.8|) it follows that J log |z — s|<i/i r (s) does 
not depend on r if z G f2oo. Indeed, the parts of /i r on the interval /fe] are the same for 
every r, and for the part v T of fi r on r r , we have by Cauchy's theorem 

/ \og{z - s)du r = — [ \og(z - s)^-ds = (]ogz)R(0) = (log*) A, (3.60) 
J 2m J Tr s 

for z G f^, which is also r-independent. The relation ( |3.60| ) is valid for any branch of 
log(z — s) which is single-valued for s in the domain enclosed by r . So J log \z — s\dfi r (s) 
is indeed independent of r, so that by ( |3.59|) , 



lim — log |P„(z)| = / log \z — s\ dfiJs) for z G f2^. (3.61) 

n^oo n J 

Since, in case L = e~ r , all zeros accumulate on T r U \fti,fa] and fi r is a probability measure 
on T r U [fa, fa], it follows from ( |3.61| ) that fi r is the asymptotic distribution of the zeros of 
P n , see [p7| , Theorem 2.3] and also p8| , Chapter III]. Note that it is important here that 
r r U [fli, fa] is the boundary of its polynomial convex hull. 

In case L — 0, it follows in a similar way that ( |1.10| ) is the asymptotic zero distribution. 



□ 
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